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ABSTRACT 


We  discuss  some  properties  of  a  certain  physically  interestiny 
nonlinear  inteyro-differential  equation  with  periodic  Doundary 
conditions.  It  is  the  natural  periodic  analoyue  of  the  intermediate 
lony  wave  equation,  and  it  provides  a  periodic  analoyue  of  the 
Benjamin-Uno  equation  in  the  appropriate  limit.  Due  to  the  speciality 
inteyral  operator,  the  equation  admits  a  Backlund  transformation,  an 
infinity  of  motion  constants,  etc.  Two  simple  periodic  solutions  are 
exhioited.  Finally  we  note  that  the  equation  .nay  be  transformed  into 
more  than  one  Kind  of  bilinear  equation. 
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The  so-called  intermediate  lony  wave  (l.L.W.)  equation  (Joseph, 


1977;  Kubota,  Ko  and  Dobbs,  1978;  can  be  written  in  the  form 


ut 


+  8uux  + 


(T  u)xx 


U 


U) 


on  -«•  <  x  <  <*>,  where 


T  u 


coth  —  (x  -  y)  u(y;  dy 
i  o 


u; 


and  the  integral  is  evaluated  in  the  principal-value  sense.  Tne 
equation  can  be  solved  on  (-«,<*>}  via  an  inverse  scatteriny  transform 
(1ST;  (Kodama,  batsuma  and  Ablowitz,  1981;,  and  has  the  analytical 
structure  associated  with  such  equations  (Joseph  and  Eyri,  ly 7o; 
batsuma,  Ablowitz  and  Koda.na,  ly79). 

The  physical  derivation  of  (1),  (8J  as  a  model  of  the  evolution 
of  long  internal  waves  of  moderate  amplitude  assumes  that  u(x;  has  a 
classical  Fourier  transform,  and  that  u  vanishes  as  |x|  +  ®.  Even  so, 
one  .nay  ask  whether  (1),  U;  admit  spatially  periodic  solutions.  This 
was  done  by  Joseph  and  Eyri  (19/8;,  Chen  and  Lee  (197y;,  and  Nakamura 
and  Matsuno  (UbU;,  using  formal  algebraic  methods.  Unfortunately, 
the  solutions  so  obtained  either  contain  errors  or  are  subject  to  a 
limitation  that  was  obscured  by  these  formal  methods. 
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Alternatively,  one  >nay  seek  an  evolution  equation  for  lony 
internal  waves  of  moderate  amplitude  that  are  spatially  periodic. 
Then  a  derivation  similar  to  the  usual  one  leads  to  (1),  but  in  the 


periodic  case  (2)  is  replaced  with 


Tu 


T(x  -  o,L)  uU)  di. 


(da; 


wnere 


^Actually,  the  physical  derivation  naturally  leads  to  the  Fourier 

representation  of  f,  given  in  (3c),  which  is  then  transformed  into 

(3b). J  Tnis  derivation  also  requires  that  J  u  dx  =  U,  which  we  nay 

-L 

impose  on  (lj  with  (3)  without  loss  of  generality.  In  (Jo),  K  denotes 
the  co.nplete  elliptic  integral  of  the  first  kind,  Z(a)  is  Jacobi's 
Zeta  function,  and  dn(a),  cs(a)  are  Jacobian  elliptic  functions. 
These  all  have  modulus  k,  determined  by  the  condition  that 
K'lk)/K(k)  =  6/L,  where  K ' ( k )  is  the  associated  elliptic  integral  of 
the  first  kind.  (All  of  these  functions  are  discussed  by  Byrd  and 
Friedman,  iy/i.)  The  purpose  of  this  note  is  to  discuss  some  of  tne 
mathematical  structure  of  (1)  with  (3;. 


An  alternative,  but  very  useful,  representation  of  T  is  its 


Fourier  series. 


T(x;  6,L)  =  i  ^  coth  p— )  exp  .  (Jc) 

n^U 

so  tnat 

Tu  =  i  coth  pY2)  un  exp  ,  V-id) 

n^U 

where  iQn}  are  the  Fourier  coefficients  of  u.  This  representation 
follows  from  the  identities  (c.f.,  Syrd  and  Friedman,  1^71,  ^yOb.ul 
and  ydd.bi; 


Z(a) 


2*  y  qm 

K  Z-j  i  .  qZm 
>T»=  1 


sin 


(t) 


cs(a)  dn(a) 


ii 

ZK 


na  £n  \  1  qin 

ZK  K  L-J  i  +  q'n 

m=i 


where  q  =  exp(-nK'/K)  =  exp(-no/l),  and  from  the  formal  representation 
(Gel'fand  and  Shilov,  1%4,  p.  JZ) 
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2 


sin  mu  . 


00 


in=l 


The  usual  operator  on  (-«,  ®)  .nay  be  recovered  simply  by  replacing  the 
sum  in  (3d)  with  an  integral  and  rescaling.  Similarly,  we  may  recover 
(2)  from  (3)  by  letting  L  -►  «,  6  fixed.  Then  k  -*•  i,  K  +  n|_/2o,  and 
one  may  show  that 


1  =,  ,  .  1  _  nx  ,  . 

—  T ( x ;  o,  L-*-)  ♦  —  coth  —  ,  (4a) 

2L  Zo  Zo 

which  reproduces  (2).  On  the  other  hand,  for  o  -*•  ®,  L  fixed,  we  nave 
k  -*■  0,  <  -*■  iT/2,  and 


f ( x ;  o+®,  l)  ♦  ‘  COt  2L  ’  Wb) 

This  is  the  well-known  Hilbert  kernel  on  (-L,L).  Wi tri  this  kernel, 
(i)  is  the  natural  periodic  extension  of  the  Benjami n-Ono  equation 
(Benjamin,  1967;  Ono,  197S).  As  one  would  expect  (1)  with  (3) 
reduces  to  the  (periodic)  KdV  equation  if  o  U,  L  fixed. 

If  u(x)  is  periodic  with  period  2L  and  with  zero  mean,  then  ( Tu j 
according  to  (2)  and  to  (3)  are  identical.  That  this  is  so  may  be 
seen  by  rewriting  (2)  in  its  Fourier  transform  representation,  and 
recalling  that  the  F.T.  of  a  periodic  function  is  a  sum  of  Dirac 
delta  functions  (Gel'fand  and  Shilov,  1964).  Thus  (1)  with  (d)  may  oe 


b 


regarded  simply  as  the  .nost  natural  way  to  write  the  l.L.w.  equation 
when  periodic  solutions  are  of  interest. 

The  operator  T  given  in  (d)  is  the  .nost  general  periodic  operator 
we  have  found  which  satisfies  the  "T-conditions" : 

T(uTv  +  vTu)  =  Tu  Tv  -  uv  ifl) 


J  (u  Tv  +  v  Tu)  dx  =  d  ,  {id) 

.00 

where  u,v  have  zero  mean.  We  find  that  (Tl),  [Id)  are  necessary 

conditions  on  the  T  operator  in  order  for  (i)  to  have  more  than  tne 
standard  number  of  conserved  quantities.  With  tnese  conditions  the 
above  evolution  equations  can  be  expected  to  be  in  the  1ST  class. 
Condition  (T2)  follows  from  the  fact  that  T(x)  is  an  odu  function.  To 
establish  (Ti)  we  use  the  representation  (dc).  Calling 

Tn  =  icoth  nuo/L,  and  un,  vn,  the  Fourier  coefficients  of  u,  v 

respectively,  assuming  ^o  =  ^o  =  ^  (i.e.,  zero  mean)  and  using  the 

convolution  theorem,  then 


b 


00  00 


T(uTv)  “  Xj  fn  Zj  ^n-m  ^n-,n  exp 


n=-»  m=-e 


Xj  «*n  £*.  f n_,n  vn-m  exp  (i~) 


,n  n 


-  H  "•«  2T  5"-» 6xp  (nr) 


n  n 


^  ^  1  t  ?  -  /i  n<tx  \ 

■  ^ (  un  l ,  vn-m  exP  (  ~J 


:n  n 


um 


r.  ^ Tn_;n  ^n"n  exp  i1^) 


,n  n 


-  u(x)  v(x)  -  r ( v Tu)  +  (TuKTv) 


where  we  have  used  the  identity 

cothA  cothii  =  1  ♦  coth^A-U)(-  cothA  +  cothli). 

The  order  of  summation  in  (b)  may  be  interchanyed  if  £| un | ,  ]T|vn| 
exist. 
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The  need  for  the  T-conditions  can  De  seen  from  the  following.  Tne 
usual  I.L.W.  equation  on(-«,»)  was  considered  by  Satsuma,  Ablowitz, 
and  Kodama  (1979).  They  showed  that  the  derivations  of  the  constants 
of  motion  and  of  the  Backlund  transformation  do  not  depend  on  the 
specific  kernel  of  T,  so  long  as  T  satisfies  conditions  (Ti)  and  (T^). 
This  implies  that  these  formulae  will  remain  valid  for  equation  U) 
with  (8)  since  the  T  conditions  are  satisfied.  Actually  one  easily 
verifies  (by  differentiation)  that  the  constants  of  motion  yiven  by 
Satsuma,  Ablowitz  and  Kodama  (1979)  are  also  constants  of  motion  of 
(1)  with  (3). 

One  may  use  the  generalized  Miura  transformation, 
u  =  ( A-t/x+iH/x+^e^)/2  to  derive  a  general  i  zation  of  the  so-called 
modified  I.L.W.  equation.  Namely  if  £[u]  =  U  represents  U),  and 
MLvj  =  t/^-  +  (1/o+a)Vx  +  (TV)XX  ^Vx  e^  +  i7x(TV)x  =  U  , 

then  using  the  T-conditions  we  have: 

Kluj  =  u MlvJ  ,  (b) 


where 


o  =  (-  ax  +  iT  dx  +  ia  e1^)  . 


Similarly  it  can  be  proven,  by  using  the  results  of  Fokas  and 
Fuchssteiner  (1980)  (whicn  also  apply  to  tne  above  ^iura-type 
transformation)  that  (i)  admits  an  auto-Backl und  transformation  if  and 
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only  if  (Tl),  (17)  are  satisfied. 

The  1ST  pair  used  by  Kodama,  Satsuma  and  ADlowitz  (l^dl)  is  also 
valid  here.  However,  the  simple  compatibility  argument  does  not 
indicate  the  importance  of  the  T-conditions.  Namely,  consider  a 
linear  scattering  problem  and  a  sequence  of  associated  time  evolutions 
of  the  form: 


iyx+  ( o  -  M  =  ^  * 


(7a; 


-  dn  v  *  n  -  l ,  Z,  . . .  1 7  b ) 

t / a j  is  to  be  thought  of  as  a  differential  Kie.nann-Hi  lbert  proDle.n  to 
find  y1  (i.e.,  y1  are  the  boundary  values  of  certain  analytic 
functions)  given  an  appropriate  font  ion  uKx).  For  each  n, 
compatibility  of  (7a, b)  yields  an  evolution  equation;  namely, 
requiring  y£t  =  y£x,  and  setting  all  coefficients  of  y+,  y",  yx,  vxx> 
etc.,  to  zero  after  using  (7a)  to  el i  innate  derivatives  of  y+,  yives 
an  algorithmic  procedure  to  determine  compatible  equations.  Tne 
equivalent  to  the  first  two  equations  of  AdV  Hierarchy  are  obtained  as 
follows.  First,  =  ax  +  A*,  whereupon  we  find  A+  -  A"  =  U. 
Taking  A1  =  A0  =  const,  the  compatible  evolution  equation  is  ut  =  ux. 
Second,  =  iax  ♦  id1^  +  iA1,  we  find  d+  -  d"  =  J, 

A+  -  A"  =  -diux.  Taking  dx  *  d0  =  i (7a.-*-1/o)  =  const, 

A1  =  iiux  -  Uo)x.  the  compatible  evolution  equation  is  (1),  without 


need  for  the  T-conditions.  The  underlying  reason  why  such  conditions 


■oust  De  added,  and  whether,  in  fact,  (d)  is  the  most  general  singular 
integral  operator  satisfying  (Tl),  (Td)  are  open  questions.  Matsuno 
( 1 yau )  has  given  a  different  algorithm  to  derive  sucn  a  hierarchy  of 
equations. 

Next  we  consider  some  special  solutions  of  (lj  with  (d).  For  a 
wave  of  permanent  forn,  at  ♦  -cax,  and  (i)  nay  be  integrated  once  to 

^  -  cj  u  +  u^  +  (Tu)x  +  A  =  0  (a) 

with  A  constant.  It  has  been  claimed  (Joseph  and  Eyri,  li>7a;  Chen 
and  Lee,  ly/yj  that 


u(x  -  ctj 


-  (nn/L;  sinh  (nno/Lj _ 

cosn(iuo/l)  +  cos  i_nn(x  -  ct;/Lj 


(yaj 


is  a  solution  of  (d).  Its  Fourier  series  representation  is 


dU) 


nn 

L 


uo 

.  ,  \  '  ,  ,  ,.n  (  mnn<^  I  mnn o  \ 

1  +  ^  2__J  i}  exp  l — cos(  ~r~) 


(ybj 


That  (J;  does  not  solve  (<3)  may  be  seen  by  computing  au/ao  from  (da) 
and  (db),  and  Tu  from  (ddj,  because  it  turns  out  that 
au/ao  =  -a(Tu)/ax.  The  correct  solution  is 


u(x  -  ct) 


_ -  (nii/L)  sinh ln»a/L) _ 

^  1 U  d  j 

cosh  (nito/L)  +  cos L(nw/Ulx  -  ct  -  iyjj 


where  y  <  -o.  Its  Fourier  series  is 


oo 

,  ,  2nn  \  .  .  ,in  /innno\  /  i  .nn»  ( <.  -  iy)\  .  . 

uU)  =  —  2_j  sinn^— —  J  exp  f - ^ - ±L\  ,  Uubj 

1 

and  its  validity  may  be  verified  by  computing  au/a6  from  u^aj  anu 
(ldbj.  This  solution  is  comp  lex- valued. 

Periodic  solutions  may  be  obtained  systematical ly  via  Hirota's 
method,  revising  slightly  the  methods  of  Cnen  and  Lee  dakamura 

and  Matsuno  (ldduj,  and  Satsuma  and  Aolowitz  Let 

t - ( x )  i  f(x+io),  let  f(x)  be  periodic  UL)  and  let  f(.zj  he  analytic  in 
the  rectangle:  -L  <  Re(z)  *  L,  -o  <■  l.n  (zj  «•  c>.  Then  by  integrating 
I  f(x-s+io)  f ( s)d»,  around  this  rectangle,  one  finds  that 


~ J  J'  T(x  -  z;  o , L j [ f + 1 z )  -  f‘(z)]dz  =  i  if+UJ  +  f‘ix)j  f  J0  , 
-L  i 


Ui; 


where  J0  is  an  unimportant  constant.  It  follows  that  if 


u(x)  =  i  l f + 1 x /  -  f“(x;j 


(U) 


where 
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<*x  (^gF1;  . 


f “(Xi  « 

and  if  f(x)  is  properly  analytic,  then  (1)  with  [d)  becomes  a  bilinear 
equation. 


+  a)  F+  • F- 


U 


Ud) 


with  the  usual  notation  Le.y.,  Dx  a.b  =  (ax-ax* )  a(xjb(x' } |  ._xj,  and 
A  =  rt(t).  We  emphasize  that  only  those  solutions  of  (id)  that  are 
analytic  in  the  rectangle  yield  periodic  solutions  of  (i),  a  fact  tnat 
wa s  overlooked  previously. 

The  simplest  real-valued  solution  of  (i;  with  (d ;  was  given  oy 
Nakamura  and  Matsuno  U^dU).  It  may  oe  written  in  the  form 


uU;  mj 


liiiili. 

L 


U- 


io);  m 


U  +  i°); 


U^aj 


where  t.  =  x-ct+xQ,  and  ,n  is  the  modulus.  The  analyticity  condition  is 
that 


±  <  < ' (m) 

L  K(m) 


U^bJ 


or,  because  o/L  =  K'(k)/Kik),  that  the  modulus  of  T  exceed  the  modulus 
of  the  solution.  There  seems  to  Oe  no  simple  formula  for  the  speed  of 
the  wave. 


Id 


Finally,  we  note  that  (13)  is  not  the  only  “bilinear"  equation 
that  may  be  obtained  from  (1),  (3).  For  example,  let  f(z)  have  a  real 
period  (2L),  and  be  analytic  in  -L  *  ReU)  *  L,  -o  *  Im(z)  <;  o  except 
for  two  poles  at  z  =  z0  and  Zq  [u  <  Im(z0)  <  oj  with  residues  b  and 
b*.  Then  the  integral  of  /  T(x-<,+  io)  f(<Jd<.  around  tne  rectanyle 
yields 


L 

f  T  ( x  -  z )  [  f + 1  z )  -  f-(z;jdz 
-  L 


i  t.f+U)  +  f"(x)  j  +  J0  - 


m 

L 


b  Tlx  - 


*o  + 


i  a)  +  o*  T  ( x  -  zj  ♦ 


■ 

io) 


(lb) 


instead  of  (11).  In  this  case,  (1^)  changes  (1)  with  (3)  into 


i  Dt  +  ~Dx  _  l>x  +  A  +  ^  LbaxT(x-z0+io)  + 


+  b< 


f(x- 


Zo  + 


i  6)  j 


F+  .  F" 


llo) 


The  analyticity  requirement  is  that  F(z)  should  be  analytic  in  the 
usual  rectangle  except  for  simple  branch  points  at  z0  and  zj$,  so  that 
f*  has  poles.  We  have  not  determined  whether  (lb;  yields  any 
solutions  of  (1),(3)  that  are  not  available  via  (13). 
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